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Abstract. - We compute the relative localization volumes of the vibrational eigenmodes 
in two-dimensional systems with a regular body but irregular boundaries under Dirichlet and 
under Neumann boundary conditions. We find that localized states are rare under Dirichlet 
boundary conditions but very common in the Neumann case. In order to explain this difference, 
we utilize the fact that under Neumann conditions the integral of the amplitudes, carried out 
over the whole system area is zero. We discuss, how this condition leads to many localized 
states in the low-frequency regime and show by numerical simulations, how the number of the 
localized states and their localization volumes vary with the boundary roughness. 



Introduction. - The problem of localization in disordered systems [1] has been subject 
to intense research for several decades. Localization of vibrational or electronical eigenmodes 
has turned out to change the physical properties, as e.g. the electric properties [1-3], the 
damping of acoustic resonators [4] and cavities [5], the band gap in semiconductors [6], the 
wave transport [7] and the density of vibrational states [8,9]. The reasons for localization 
have been discussed for long and it is generally accepted that localization always involves 
some geometrical or structural irregularity of the system. This irregularity can be present 
as bulk irregularity or as boundary roughness. The first case has been widely investigated 
and standard models are e.g. the Anderson model [1,3] and the percolation model [10]. 
Less works however addressed the problem of localization in systems with an ordered bulk 
material and an irregular boundary. One model for this problem is the model of fractal drums, 
where localized states have been found under Dirichlet [9, 12] as well as under Neumann [8,9] 
boundary conditions. 

Recently [11], localized states have also been found in 2d-systems with non- fractal bound- 
ary roughness under Neumann boundary conditions, while under Dirichlet conditions they are 
rare and seem to be linked to the existence of confined regions (which indeed occur in the case 
of fractal drums [12] or in systems with hard scatterers [13]). In this Letter we focus on this 
phenomenon and consider several types of systems with non- fractal irregular surfaces. Under 
Neumann boundaries we find many localized states in all systems with arbitrarily shaped 
irregular boundaries. We elucidate this behavior by numerical simulations and show how it 
can be explained by using a sum-rule for the Neumann case. 
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Fig. 1 - Localized modes of several systems with rough boundaries under Neumann boundary condi- 
tions for different system geometries, (a-b) L-shaped system, (c-d) triangular systems and (e-f) tooth 
systems. The amplitudes are indicated by different gray levels. The white and black regions stand for 
positive and negative amplitudes, respectively and the neutral gray tone stands for nearly zero ampli- 
tude. The eigenfrequencies and relative localization volumes are (a) w » 0.060 (fc/m) 1/2 , Viae ~ 0.13, 
(b) w ~ 0.088 (k/m) 1/2 , V ioc « 0.13, (c) u w 0.125 (fc/m) 1/2 , Vi oc » 0.07, (d) u> « 0.226 (fc/m) 1/2 , 
Vioc » 0.06, (e) u « 0.072 (fc/m) 1/2 , Hoc ~ 0.17 and (f) uj 2 » 0.083 (fc/m) 1/2 , Vice ~ 0.09. 



Localized Modes. - We start with the Hclmholtz equation for the vibrational amplitude 
ip a (x, y) of a membrane that is located in the xy-plane and vibrates in z-direction, 

uj 2 

Aip a (x,y) = i> a (x,y), (1) 

with the state index a, the eigenfrequency 07 Q and the sound velocity c. The boundary may 
have an arbitrary shape and may either be fixed (Dirichlet case) or vibrate freely (Neumann 
case). The latter case applies e.g. for gas vibrations in acoustic resonators, where ip describes 
the pressure distribution of the gas [5] . For calculating the eigenfunctions and -frequencies, we 
discretize Eq. Q on a square lattice, where nearest-neighbor masses m are coupled by linear 
springs k and diagonalize it by the Lanczos algorithm. When the number of masses tends to 
infinity, we approach a continuous system, where a 2 k/m — * c 2 , with the lattice constant a. 

For regular systems, all states are extended for both types of boundary conditions. The 
eigenstates of a rectangle of side lengths L x and L y , e.g., are simple sine- or cosine functions 
under Dirichlet and under Neumann boundary conditions, respectively. For a discretized 
rectangle of Z x = L x /a and Z y = L y /a lattice segments along the x- and the y-axis we have 
w n,m = (2k/m)(2- cos(nw/£ x ) - cos(rmr/£ y )) where li = Zi + 1, n £ [0,Z X ], m £ [0,Z y ] 
under Neumann and £j = Zi, n £ [1,Z X — 1], m £ [l,Z y — 1] under Dirichlet boundary 
conditions [8,14]. In the continuous limit n <C Z x ,m <C Z y the eigenvalues converge to 
< m = vr 2 c> 2 /^ + m 2 /^). 

For systems with irregular boundaries, the eigenstates show an irregular pattern of several 
mountains and valleys of similar but not identical shapes, whose extensions refer roughly to 
half a wavelength A. Therefore, at higher frequencies the mountains and valleys are narrower 
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Fig. 2 - Relative localization volumes Vioc versus lu 2 for the systems of Fig. under (a) Dirichlet and 
(b) Neumann boundary conditions. The symbols o, □ and o stand for the geometries of Fig. H^a-b) 
(L-shaped system), (c-d) (triangular system) and (e-f) (tooth system), respectively. The straight line 
refers to the value of 4/9 of an ordered square membrane. 



and therefore more numerous than at low frequencies. This can be seen in Fig. QJa-f), where 
some representative localized states of different systems are shown. For the localized states 
of Fig. 1, the absolute values of the vibrational amplitudes are only large at one part of the 
system and decay rapidly towards the other system side. 

In order to quantify the degree of localization, we determined numerically the relative 
localization volumes VJ " (participation ratio) [15], 

v (a) = YsL= I 

lo ° A A J \tp( a )\ 4 dxdy 

where ip is normalized according to J \ip(x, y)\ 2 dxdy = 1 and A is the membrane area. V\ oc is 
equal to one for constant functions, whereas for regular sine or cosine-functions of an ordered 
square membrane we have V\ oc = 1 for n = m = 0, V\ oc = 2/3 for n or m — and V\ oc = 4/9 
for n,m > 0. For a Gaussian function, we have V\ oc = 1/3. If Vioc is much smaller than 
this value, we call the function localized. For comparison, the values of V\ oc are given in the 
caption of Fig. ^ In Fig. [3 we show all values of V\ oc of the low-frequency modes of the 
systems of Fig. \T\ior (a) Dirichlet and (b) Neumann boundary conditions. In Fig.^a) small 
values Vioc below 0.2 do not occur and most values are above 0.33. In Fig.[5{b) on the other 
hand, the distribution of the localization volumes is much broader and shows a tail, with a 
large number of values even below 0.1. So, there must be a strong localization mechanism in 
the Neumann case, which apparently does not exist under Dirichlet boundary conditions. 

To understand this difference, we recall that a sum rule applies in the Neumann case for 
each wavefunction (except for the first one, where ui a = 0) [8] . This can be seen by integrating 
Eq. ||TJ over dx dy and applying the divergence theorem, yielding — (w^/c 2 ) J J dxdy ip a (x, y) = 
J r (Vip a )nds, where T is the boundary of the membrane, n its normal component and the 
double integral on the left-hand side is carried out over the membrane area. With the Neumann 
condition dip/dn — (Vip)ri — along T, we find the sum rule 

dxdy ip(x, y) — for w^0, (3) 

which is rather restrictive. As we will see in the following it strongly limits the way, how 
a wavefunction that is well-suited to a given rough boundary evolves into the rest of the 
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Fig. 3 - The L-shaped system AB. (a) The geometry, (b) 20th Eigenmode of the total system AB 
in part A alone, (c) Eigenmode (n, m) = (0,2) of the isolated rectangle A. (In (b,c), part A is 
shown enlarged as compared to (a).) (d) Example of the 6th eigenmode under Dirichlet boundary 
conditions, arising from the (1, l)-mode of rectangle A for comparison. 

system and is the key in understanding the occurence of the localized states under Neumann 
conditions. In the Dirichlet case, where Eq. does not apply, the absolute value of the sum 
/ / dxdyip(x,y) can take any value, i.e. the mountains and valleys of the mode do not need 
to cancel each other. 

The L-shaped system. - To show, how the sum rule @ leads to localized states, we first 
consider the L-shaped system, which is composed of two rectangles A and B of sizes A x x Ay 
and B x x By, respectively (see Fig. EH a)). We assume that A x and A y are no divisors of B x 
and By (see below), i.e. B x ^ iA x and B y ^ jA y , i,j € N. Figure|3fb) shows the eigenstate 
of Fig. [IJa), in A alone (zoomed to a larger size). Comparing it to the state (n,m) = (0,2) 
of the isolated rectangle A (cf. Fig. 02c)), we see that apart from some decay from the left 
towards the right side of the system in (b), both mode patterns are quite similar. Also the 
frequencies u> of both states are very close, as can be seen in Table 1 ((n, m) = (0, 2), a = 20). 
Also many other rectangular eigenstates, with n = 0, m > show up in the spectrum of the 
L-shaped system with nearly zero amplitudes in region B [16]. 

We now show how the localization is related to the sum rule J2Jl: If modes, originating 
from regular modes in A are extended, they must evolve into region B. An integer number of 
the mountains and valleys, that are optimized for region A must travel to region B, satisfying 
(i) Eq. at each lattice point and (ii) the sum rule 10 ■ This can be achieved as follows: 
The integral Sa = JaJ dxdy tp(x, y) over region A alone is close to zero, as it arises from an 
eigenmode of rectangle A. In region B, the corresponding integral Sb = J B J dxdy ip(x,y) 
must therefore just compensate the value of Sa, i-e. Sb = —Sa- As Sa is small, it is clear 
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Table I - Comparison of the eigenfrequencies uj„ im (in units of (k/rn) 1 ^ 2 ) of rectangle A with the side 
lengths (numbers of particles) N x = 79, N y = 110 and ui a of the L-shaped system for (i) Neumann 
and (ii) Dirichlet boundary conditions. The last columns show Vk>c,a of the L-shaped system. 
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that this condition is easier to fulfill, when the absolute values of ip in region B are small. 
When the \tp(x, y)\ are large in region B, the single mountains and valleys need to cancel each 
other accurately, which is in most cases impossible because of the mismatch between the two 
rectangles A and B. If on the other hand the eigenfunction decays to very small values in 
region B, also the integral Sb stays small. Moreover, by a suited decay rate it can always 
be adjusted such that Sb = — Sa, thereby fulfilling the sum rule ©. This clearly leads to 
localized modes. 

Under Dirichlet boundary conditions, the situation is different, as the sum rule (|3Jl is not 
valid and therefore the possible modes are much less restricted. One typical mode is shown 
in Fig. Old), which arises from the (n,m) = (1, l)-mode of rectangle A. Obviously, this mode 
is extended and clearly, J J ip(x,y)dxdy is unequal to zero. 

To test this assumption, we analyzed the mode pattern and the eigenvalues of the modes, 
arising from rectangle A under (i) Neumann and (ii) Dirichlet conditions. In Table 1, we 
compare their frequencies u a and relative localization volumes V\ oc a to the corresponding 
frequencies u> n , m of rectangle A. In the Neumann case, most of them are indeed localized. 
(Not shown here, there are also localized modes of slightly higher localization volumes arising 
from the other two possible rectangles of the upper and the lower part of the system.) Most 
Dirichlet modes, on the other hand, possess much higher values of Voc- 

Localization in general systems. - We now consider more complex systems, as e.g. the 
systems of Fig.^c-f). The localization in these cases can be understood by the same mecha- 
nism as above, i.e. by eigenmodes of some boundary subsystems whose amplitudes decrease 
when entering other regions of the system in order to fulfill the sum rule (JSJ. This can be seen 
in Fig. ^ where we can most easily identify rectangular eigenfunctions in boundary structures 
(c-d), or modes of several close-lying boundary rectangles that couple (e-f). 

We now discuss, how the number of localized states and the values of the smallest localiza- 
tion volumes can be influenced by the boundary. To this end, we calculated the histograms of 
the Moc of several systems under both types of boundary conditions in the frequency regime 
of < ijj 2 < QAk/m. Figure a) shows as an example the complete histogram of the tooth 
system under Neumann conditions, which has a strong maximum at V\ oc ~ 0.33, which refers 
to a Gaussian distribution of the amplitudes and does not change when going to other systems 
or boundary conditions. 

We are mostly interested in the tails of the histograms towards smaller values, which 
depend sensibly on the boundary details. In Fig. 0fb-e) we show the histogram tails for 
< Moc < 0.3 for (b) the L-shaped system, (c) the triangular system, (d) the tooth system 
and (e) a "confined" L-shaped system where region A was made much longer and narrower. 
The histograms for the Neumann states are indicated by white columns and the histograms 
for Dirichlet states by black columns. 

We can see that in the first three cases, localized states under Neumann conditions are 
much more frequent than under Dirichlet conditions. In the normal L-shaped system we 
found only two modes with Vi oc < 0.2 in the given frequency range in the Dirichlet case, 
but a large number of localized states in the Neumann case with Vioc going down to values 
even smaller than 0.1. This phenomenon is even stronger in the case of the triangular system 
of Fig. life) and most pronounced in Fig. Efd), which corresponds to the system with many 
narrow boundary "teeth" and where no localized states are found in the Dirichlet case. 

In the Neumann case, the number and Voc-values of the localized states change with the 
shape of the systems. First, as complex boundaries possess more small substructures, where 
localized eigenfunctions can exist, the number of localized states increases with the complexity 
of the boundary, i.e., the systems of Figs.Qfc) and (d) possess more localized states than the 
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Fig. 4 - (a) Complete histogram of the tooth system under Neumann boundary conditions in the 
frequency regime < lj 2 < OAk/m. One can recognize the maximum of the histogram close to the 
modes with Gaussian distributions, i.e. with Vioc ~ 0.33. (b-e) Tails of the histograms of Vioc in the 
range between V\ c = and 0.3 for (b) the L-shaped system, (c) the triangular system, (d) the tooth 
system and (e) the confined L-shaped system with a much longer and narrower A region. White 
columns refer to Neumann, black columns to Dirichlet boundary conditions. 



L-shaped systems. (Note the different range of the y-axis.) Second, the smallest values of 
Vioc decrease with the relative size of the boundary substructures, i.e. the values of V\ oc in 
Figs. EIc,d) decrease to smaller values than in (b,e). Third, also the total system area A 
influences the relative number of localized states. Whereas the number of extended states 
grows cx A, the number of substructures and thus of the localized states grows cx \f~A~ (given 
a constant density of substructures along the boundary). Therefore, the relative number of 
localized states should become smaller with increasing system size. Nevertheless, the number 
of localized states should stay non-negligiable in most cases of rough or porous materials that 
are built up by smaller clusters that do not grow to arbitrarily large sizes. These relations in 
mind, one could use the localization behavior also for technical applications, as e.g. acoustic 
damping in systems where the surface is covered with an absorbing material. It has already 
been shown that localized states are much more strongly damped than extended states [4,5], 
so that a suitable choice of the size of the substructures could be very efficient in damping 
the modes of a certain frequency range. 

Only the confined L-shape system of Fig. 0Je) is exceptional, as it possesses localized 
states with values of V\ oc close to 0.1 also under Dirichlet conditions. This is in agreement 
with earlier numerical and experimental results on fractal drums [12] and on systems with 
hard scatterers [13], where localized states under Dirichlet conditions have also been found 
in confined regions. Clearly, localization under Dirichlet conditions is not forbidden, but it is 
quite rare and the localization mechanism is not the same as the one described here. Whereas 
under Neumann conditions, slight boundary irregularities are sufficient to create localized 
states via the sum rule 0, under Dirichlet conditions, localization seems to be linked to the 
occurence of confined regions in the system. 
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Summary and Conclusions. - In summary, we have investigated systems with non-fractal 
rough boundaries and computed the relative localization volumes V\ oc of the eigenstates. We 
found that under Neumann boundary conditions, many localized states are coexisting with 
extended states in the low-frequency regime, whereas under Dirichlet conditions localized 
states are very rare. This difference can be explained by a sum rule that only applies under 
Neumann boundary conditions and implies a localization mechanism for systems with rough 
boundaries. By investigating systems of different shapes, we found that number and relative 
localization volumes of the localized states can be triggered by the system geometry, i.e. by 
the complexity of the boundaries and the size ratio of the body and the substructures. 

The above arguments do not apply, when the side lengths of the boundary rectangles are 
integer multiples of the main body and the mountains and valleys of the boundary modes 
fit also into the body and therefore cancel under the integral ©, even without any decay. 
However, in natural systems, where the boundary structures are no simple rectangles but 
random structures, this case will hardly occur. 

The sum rule is easier to fulfill in the high-frequency range, where the mountains 
and valleys are small and numerous. Accordingly, the described localization mechanism only 
applies below some limiting frequency regime whereas in the high-frequency regime, the local- 
ization volumes are close to the Gaussian value of Vi oc « 0.33. There is no crossover between 
the regimes, but rather a disappearance of the localized ones around frequencies between 
1.5k/m < Lo 2 < 2k/m, which refer to a half-wavelength A/2 w 2 — 3 lattice constants a, the 
exact value depending on details of the system. This dependence on a in discretized systems is 
straightforward, as a is the smallest common divisor of the different boundary substructures 
and the main body of the system. In continuous systems, we expect that the localization 
mechanism also terminates at wavelengths comparable to the value of the smallest common 
divisor of the different lengths of the boundary roughness. 
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